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Abstract 

We analyze combinatorial optimization problems over a pair of random point sets {X, y) in R'* of 
equal cardinal. Typical examples include the matching of minimal length, the traveling salesperson 
tour constrained to alternate between points of each set, or the connected bipartite r-regular graph 
of minimal length. As the cardinal of the sets goes to infinity, we investigate the convergence of such 
■ bipartite functionals. Q 

(N 

O ■ 1 Introduction 

o 

This work pertains to the probabilistic study of Euclidean combinatorial optimization problems. The 
starting point in this field is the celebrated theorem of Beardwood, Halton and Hammersley [2] about 
the traveling salesperson problem. Its ensures that given a sequence {Xi)i>i of independent random 
variables on M."^, d > 2 with common law jj, of bounded support, then almost surely 

CM ■ 



in 



X 



lim n^-iT(Xi,...,X„) =/3rf / /^'i 

n— >oo J 



Here (3d is a constant depending only on the dimension, / is the density of the absolutely continuous 
part of /.i and 

n-l 

<N : T(Xi, ...,Xn)= in| ^ - + - X,(„)| 

^■f-^ ■ is the length (for the canonical Euclidean distance) of the shorstest tour through the points Xi, . . . , X„. 

I In the above formula 5„ stands for the set of permutations of {1, 2, . . . , n}. Very informally, this result 

O^l ■ supports the following interpretation: when the number of points n is large, for /x almost every x, if 

I the salesperson is at V^ = a; then the distance to the next point in the optimal tour is comparable to 

■ (3{d){nf{x))^^^'^ if f{x) > and of lower order otherwise. This should be compared to the fact that the 
distance from V^ = a; to {Xj,j < n and j ^ i} also stabilizes at the same rate. 

Later, Papadimitriou [5] and Steele [T^ have initiated a general theory of Euclidean functionals 
i^({Vi, . . . ,Xn}) that satisfy almost sure limits of this type. We refer the reader to the monographs of 
Steele [l5j and Yukich ^9] for a full treatment of this now mature theory, and present a short outline. 
It is convenient to consider multisets rather than sets, so throughout the paper {xi, . . . will stand 

■ for a multiset (the elements are unordered but may be repeated). The umbrella theorem in |19j puts 
forward the following three features of a functional F on finite multisets of R"^: 

• _F is 1-homogcneous if it is translation invariant and dilation covariant: 

F{a + XX) = XF{X) 

for all finite multisets X, all a € M'* and A G M+. 

• The key assumption is subadditivity: F is subadditive if there exists a constant C > such that 
for all multisets X,y in the unit cube [0, 1]^, 

F{x uy)< F{x) + F{y) + c. 



AMS 2000 Mathematics Subject Classification : 90C27 (Primary); 39B62 ; 60D05 (Secondary). Keywords : combi- 
natorial optimization, minimal matching, geometric probability. 



1 



As noted by Rhcc in [I^], this assumption implies that there is another constant C" such that for 
ah multiset in [0, l]'^, 

\F{X)\<C'{caTd{X))^-i . (1) 

Hence the worst case for n points is at most in n^~^ and the above mentioned theorems show that 
the average case is of the same order. 

• The third important property is smoothness (or regularity). A functional F on finite multisets M** 
is smooth if there is a constant C" such that for all multisets X,y,Z in [0, 1]'', it holds 

\F{X uy)- F{X U Z)\ < C" (card(3^)i~3 + cai-d{Z)^'^^ . 

These three properties are enough to show upper limits for F, on the model of the Beardwood, Halton, 
Hammersley theorem. To have the full limits, the umbrella theorem of |19| also requires to check a few 
more properties of a so-called boundary functional associated with F. 

Next, let us present a classical optimization problem which does not enter the above picture. Given 
two multi-subsets of R'^ with the same cardinality, X = {Xi, . . . , X„} and y = {Yi, . . . , F„}, the cost of 
the minimal bipartite matching of X and y is defined as 

n 

M,{X,y) = min V|X,-F,(,)|, 

CTGOn 

4=1 

where the minimum runs over all permutations of {1, . . . , n}. It is well-known that n~^Mi {\Xi\^^i, {Yi}"=i) 
coincides with the power of the Li-Wasserstein distance between the empirical distributions 

\ n ^-^ n ^-^ I 

i i 

hence it is easily seen to tend to 0, for example when has bounded support. Recall that given two 
finite measures /ii on with the same total mass, 

W^i(^i,/i2)= inf / \x - y\d-K(x,y), 

where n(/^i, is the set of measures on (R'')2 having as first marginal and [I2 as second marginal 
(see e.g. [TU1[TH] for more background). Note that for all finite multisets A", y in [0, 1]'' with card(A') = 
card(3^), 

Mx{X,y) < Vrfcard(A'), 

and equality holds for some well-chosen configurations of any cardinal (all elements in A" at (0, • • • , 0) and 
all elements in 3^ at (1, • • • , 1)). Hence, an interesting feature of L (as well as others bipartite Euclidean 
optimization functionals) is that the growth bound assumption ([1]) fails, hence it is not subadditive in 
the above sense. However Dobric and Yukich have stated the following theorem: 

Theorem 1 ([4]). Let d> 3 be an integer. Assume that fi is a probability measure on R"^ having a bounded 
support. Consider mutually independent random variables {Xi)i>i and {Yj)j>i having distribution ^. 
Then, almost surely, 

limn^-iMi({Xi,...,X„},{Yi,...,y„}) =/?i(d) / f^, 

where f{x)dx is the absolutely continuous part of fi and I3i{d) is a constant depending only on the 
dimension d. 

When / is not the uniform measure on the unit cube, there is an issue in the proof of [4] that 
apparently cannot be easily fixed (the problem lies in their Lemma 4.2 which is used for proving that 
the liminf is at least /3i((i) Jj^^ f~^)- In any case, the proof of Dobric and Yukich is very specific to 
the bipartite matching as it uses from the start the Kantorovich-Rubinstein dual representation of the 
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optimal transportation cost. It is not adapted to a general treatment of bipartite functionals. The 
starting point of our work was recent paper of Boutet de Monvel and Martin [5] which (independently of 
[3]) establishes the convergence of the bipartite matching for uniform variables on the unit cube, without 
using the dual formulation of the transportation cost. Building on their approach we are able to propose 
a soft approach of bipartite functionals, based on appropriate notions of subadditivity and regularity. 
These properties allow to establish upper estimates on upper limits. In order to deal with lower limits we 
adapt to the bipartite setting the ideas of boundary functionals exposed in [TH] . We are able to explicitly 
construct such functionals for a class of optimization problems involving families of graphs with good 
properties, and to establish full convergence for absolutely continuous laws. Finally wc introduce a new 
notion of inverse subadditivity which allows to deal with singular parts. 

This viewpoint sheds a new light on the result of Dobric and Yukich, that we extend in other respects, 
by considering power distance costs, and unbounded random variables satisfying certain tail assumptions. 
Note that in the classical theory of Euclidean functionals, the analogous question for unbounded random 
variables was answered in Rhee [T3| and generalized in [T^ . 

Let us illustrate our results in the case of the bipartite matching with power distance cost : given 
p > and two multi-subsets of R'', X ~ {^i, . . . , Xn} and y = {Yi, . . . , y„}, define 



Mp{X,y) = min^|X,-y,(,)|P, 



=1 

where the minimum runs over all permutations of {1, . . . , n}. Note that wc have the same result for 
the bipartite travelling salesperson problem, and that our generic approach puts forward key properties 
that allow to establish similar facts for other functionals. As mentioned in the title, our results apply 
to relatively high dimension. More precisely, if the length of edges are counted to a power p, our study 
applies to dimensions d > 2p only. 

Theorem 2. Let < 2p < d. Let fi he a probability measure on M.'^ with absolutely continuous part 
f{x) dx. We assume that for some a > -j^i^, 

\x\"d^{x) < +00. 

Consider mutually independent random variables {Xi)i>i and (^)j>i having distribution fi. Then there 
are positive constants j3p{d), I3p{d) depending only on {p,d) such that the following convergence holds 
almost surely 

limsupn5-iAfj,({Xi,...,X„},{yi,...,r„}) < /3p(d) / 

n jRd 

liminf7i5-iA/p({Xi,...,X„},{yi,...,r„}) > 13' (d) [ /i-^ 

Jr" 

Moreover 

limnS-iMp({Xi,...,X„},{ri,...,y„}) / f'-^ 

provided one of the following hypothesis is verified: 



p is the uniform distribution over a bounded set Q C M. with positive Lebesgue measure 



• d € {1,2}, p £ (0, (i/2) or d > 3, p £ (0, 1], and f is up to a multiplicative constant the indicator 
function over a bounded set fl C M.'^ with positive Lebesgue measure. 

Our constant (d) has an explicit expression in terms of the cost of an optimal boundary matching 
for the uniform measure on [0, 1]'' (see Lemma [57]). We strongly suspect that f3p{d) = f3p{d) but we have 

not been able to solve this important issue. Also, assuming only a > J^^^ , we can establish convergence 
in probability. A basic concentration inequality implies that if fi has bounded support the convergence 
holds also in L'' for all q > 1. 

The paper is organized as follows: Section [2] presents the key properties for bipartite functionals (ho- 
mogeneity, subadditivity and regularity) and gathers useful preliminary statements. Section[3]establishes 
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the convergence for uniform samples on the cube. Section S] proves upper bounds on the upper Hmits. 
These two sections essentially rely on classical subadditive methods, nevertheless a careful analysis is 
needed to control the differences of cardinalities of the two samples in small domains. In Section [SJ we 
introduce some examples of bipartite functionals. The lower limits are harder to prove and require a 
new notion of penalized boundary functionals. It is however difficult to build an abstract theory there, 
so in Section |51 we will first present the proof for bipartite matchings with power distance cost, and 
put forward a few lemmas which will be useful for other functionals. We then check that for a natural 
family of Euclidean combinatorial optimization functionals defined in i j5.3l the lower limit also holds. 
This family includes the bipartite traveling salesman tour. Finally, Section [7] mentions possible variants 
and extensions. 

2 A general setting 

Let Aid be the set of all finite multisets contained in W^. We consider a bipartite functional: 

L : Md X Md ^ K+. 

Let p > 0. We shall say that L is p-homogencous if for all multisets X,y, all a e R'' and all A > 0, 

L{a + XX,a + Xy) = XPL{X,y). (Hp) 

Here a + X{xi, . . . ,Xk} is by definition {a + A.ti, . . . , a + Xxk}. For shortness, we call the above property 
CHp). Note that a direct consequence is that i(0,0) = 0. 

The functional L satisfies the regularity property (TZp) if there exists a number C such that for all 
multisets X ,y, Xi,yi, X2,y2, denoting by A the diameter of their union, the following inequality holds 

L{X UXi,yu yi) < L{X UX2,yu 3^2) + CAP(card(A'i) + card(A'2) + card(3^i) + card(3^2)) • (7^p) 

The above inequality implies in particular an easy size bound: L{X, y) < CA^'(card(A') + card(3^)) when 
L(0,0) = 0. 

Eventually, L verifies the subbaditivity property (Sp) if there exists a number C such that for every 
k > 2 and all multisets {Xi,yi)f^i, denoting by A the diameter of their union, the following inequality 
holds 

i( U U ^ I] W^y^) + + IcardCA-,) - card(J^,)|) • (Sp) 

i—1 i—1 i—1 i—1 

Remark 1 . A less demanding notion of " geometric subadditivity" could be introduced by requiring the 
above inequality only when the multisets Xi U yi lie in disjoint parallelepipeds (see [19] where such a 
notion is used in order to encompass more complicated single sample functionals). It is clear from the 
proofs that some of our results hold assuming only geometric subadditivity (upper limit for bounded 
absolutely continuous laws for example). We will not push this idea further in this paper. 

We will see later on that suitable extensions of the bipartite matching, of the bipartite traveling 
salesperson problem, and of the minimal bipartite spanning tree with bounded maximal degree satisfy 
all these properties. Our main generic result on bipartite functionals is the following. 

Theorem 3. Let d > 2p > and let L be a bipartite functional on with the properties (Hp), (T^p) 
and {Sp). Consider a probability measure ^ on R'* such that there exists a > f^2p ^^^'^ 

J \x\"dp.{x) < +00. 

Consider mutually independent random variables {Xi)i>i and {Yj)j>i having distribution fi. Let f be a 
density function for the absolutely continuous part of fi, then, almost surely, 

L{{Xi,- ■ ■ ,Xn},{Yi,- ■ ■ ,Yn}) f rl-E, 

hmsup — ^ < /3l / / ^ 

n->-oo n d J 
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for some constant /3l depending only on L. Moreover, if fi is the uniform distribution over a bounded 
set il with positive Lebesgue measure, then there is equality: almost surely, 

n->oo fi} d 

Beyond uniform distributions, lower limits are harder to obtain. In Sectional we will state a matching 
lower bound for a subclass of bipartite functionals which satisfy the properties {Hp), {Tip) and {Sp) (see 
the forthcoming Theorem 1361 and, for the bipartite traveling salesperson tour, Theorem I37p . 

Remark 2. Let 5(1/2) = {x g R'' : |x| < 1/2} be the Euclidean ball of radius 1/2 centered at the origin. 
It is immediate that the functional L satisfies the regularity property {TZp) if it satisfies property {'Hp) 
and if for aU multisets X,y,Xi,yi,X2,y2 hi B{l/2), 

L(X \JXi,y\J yi) < L{X \JX2,y\J 3^2) + C(card(<Yi) + card(A'2) + card(3;i) + card(3;2)). (7^) 

Similarly, L will enjoy the subbaditivity property {Sp) if it satisfies property {Hp) and if for every fc > 2 
and aU multisets (A;,3^,)JLi in 5(1/2), 

L(\jx,,{jy,^ <5I^('^»'3^0 + C^I](l+|card(A',)-card(3;,)|)- (5) 

2—1 2—1 i—1 i—1 

The set of assumptions {Hp), {Tip), {Sp) is thus equivalent to the set of assumptions {Hp), (72,), {S). 

2.1 Consequences of regularity 
2.1.1 Poissonization 

For technical reasons, it is convenient to consider the poissonized version of the above problem. Let 
(Xj)j>i, {Yi)i>i be mutually independent variables with distribution ^. Considering independent vari- 
ables A^i, N2 with Poisson distribution 'P{n), the randoms sets {Xi, . . . ,Xn^} and {Yi, . . . ,¥^2} are 
independent Poisson point processes with intensity measures n^. For shortness, we set 

L{n^l) L{{X,, ...,Xn,}. {Yi, ■ • ■ , ^Vj). 

When dfj.{x) = f{x)dx we write L{nf) instead of L{n^). Note that whenever we are dealing with 
Poisson processes, n € (0,-|-oo) is not necessarily an integer. More generally L(i^) makes sense for any 
finite measure, as the value of the functional L for two independent Poisson point processes with intensity 
v. 

Assume for a moment that the measure /x has a bounded support, of diameter A. The regularity 
property ensures that 

\L{{X,, . . . {n, . . . , r„}) - L({Xi, . . . ,XArJ, {Fi, . . . ,yArJ)| < CAP{\N, - n| + {N^ - n\). 

Note that ElTVi-nl < {E{Ni-n)^)^^^ = Var(7V,) = y/Ti. Hence the difference between EL ({X,}f^ 1, {yjf^J 
and KL{nii) is at most a constant times ^/n = o{n^~P/'^) when d > 2p. Hence in this case, the orig- 
inal quantity and the poissonized version are the same in average at the relevent scale n^~P/'^. The 
boundedness assumption can actually be relaxed. To show this, we need a lemma. 

Lemma 4. Let a > 0, n > and let ^ be a probability measure on K*^ such that for all t > 0, 
^jl{{x; \x\ > t}) < ct^°'. Let X, y be two independent Poisson point processes of intensity and 
Tn = max{|.Z| : Z Cz X U y}. Then, for all < j < a there exists a constant K = K{c, a, 7) such that 
for all n > 1, 

E[rj]T < Kni. 

Moreover the same conclusion holds if X ~ {Xi, . . . , Xn}, y ~ {Yi, . . . , Yn} are two mutually indepen- 
dent sequences of n variables with distribution /z. 
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Proof. For t > 0, let At = {x E R'^ : \x\ > t} and g{t) = J^^ dfi. By assumption, fi{At) < ct'". We start 
with the Poisson case. Since X, y are independent, we have P(T„ < t) = P{X (1 At = 0)^ = e-^^^C^t), 
Therefore, using 1 — e~" < min(l,M), 



poo 

E[rj] 7 / t'^-ip(r„ > t)dt 

Jo 



Jo 

a — 7 

For the second case, since P(T„ > t) = 1 — (1 — < min(l, 2n/i(At)) the same conclusion holds. □ 

Proposition 5. Lei d > 2j3 > 0. Let /i fee a probability measure on such that J \x\" d^{x) < +oo for 
some a > ^^2p ■ {Xi)i>i, (li)i>i fee mutually independent variables with distribution /i. If L satisfies 
the regularity property (TZp) then 

j.^ m{x^n=l,{y^}?=l)-mnf^) ^ ^ 

Remark 3. Wc have not yet proved the finitcness of the above integrals. This will be done later. The 
proof below show that the expectations are finite at the same time. So the above statement is established 
with the convention oo — oo = 0. 

Proof. Let A'^i and A^2 be Poisson random variables with mean value n. Let T = max{|Z| : Z G 
{Xi,--- ,Xn,} U {n,--- ,Yn,}} and S = max{|Z| : Z € {Xi,--- ,X„} U {Fi,--- ,y„}}, with the 
convention that the maximum over an empty set is 0. The regularity property ensures that 

\L{{Xi, . . . {Fi, . . . ,y„}) - ...,XnA, {Yi, . . . , FjvJ) I < CiT + Sf (|7Vi - n\ + \N2 - n\) . 

Taking expectation gives, using Cauchy-Schwarz inequality and the bound (a+fe)'' < max(f , 2i^^){a'^+b'^) 
valid for a, fe, g > 

|EL({Xi, . . . ,X„}, {n, . . . , r„}) - L({Xi, . . .,Xn,}, {Yi, . . . , YatJ) I 



E[T^P] + E[S'^P]j ' (e[\Ni - ] + E[\N2 - n\ 
2^(e[T^p] + E[S^P] 



Since a > 2p, by LemmaH for some c> and ah n > 1, E[T'^p] < cn'^P/" and E[S^p] < cn^Pl°'. Hence 
the above difference of expectations is at most a constant times na"*""? , which is negligeable with respect 
to v}~~i since a is assumed to be large enough. □ 



2.1.2 Approximations 

Proposition 6. Assume that a bipartite functional L satisfies the regularity property (TZp). Let m,n > 
and jjL be a probability measure with support included in a set Q. Then 

EL(n/i) < EL{m^) + Cdiam(g)P|m - n\. 

Proof. Assume n < m (the other case is treated in the same way). Let (Ai)i>i, (yi)i>i, A''i, iV2, i^i, 7^2 
be mutually independent random variables, such that for all i > 1, Aj and Yi have law ^, and for 
i £ {1, 2}, the law of Nj is P{n) and the law of Kj is Vim — n). Then Mi = Ni + Ki is P(m)-distributed. 
Then {Ai, . . . , A^Tj} and {Yi, . . . ,^^3} are independent Poisson point processes of intensity nfi, while 
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Then {Xi, . . . , and {Yi, . . . , YM2} are independent Poisson point processes of intensity m/i. By 

the regularity property, 

L{{Xi, ...,Xn,}, {Yi, . . . , FjvJ) < L{{Xi, . . .,Xn,+k,}, {Yi,-- ■ ,>V.+kJ) + Cdiam(Q)P(A-i + K2). 

Taking expectations gives the claim. □ 

Applying the above inequality for m = gives a weak size bound on EL(i/). 

Corollary 7. Assume that L satisfies (TZp) and L{^,^) ^ Q (a consequence of e.g. (Hp)), then if v is a 
finite measure with support included in a set Q, 

EL{v) < Cdiam(Qf i/(g). 

Recall the total variation distance of two probability measures on R"* is defined as 

dTv{fi,fJ.') = sup{\n{A) - fi'{A)\ : A Borel set of R''}. 

Proposition 8. Assume that L satisfies (TZp). Let be two probability measures on R"* with bounded 
supports. Set A be the diameter of the union of their supports. Then 

EL{n^i) < EL{nfi') + iCAP ndTviii, n'). 

Proof. The difference of expectations is estimated thanks to a proper coupling argument. Let tt be 
a probability measure on R'^ x R'' having /i as its first marginal and as its second marginal. We 
consider mutually independent random variables Ni, N2, {Xi, X-)i>i, {Yi,Y-)i>i such that Ni,N2 are 
V{n) distributed and for all i > 1, {Xi,Xl) and (Yi,Y-) are distributed according to tt. Then the 
random multisets 

X ^{Xi,...,XnA and 3^ = {Fi, . . . , FatJ 

are independent Poisson point processes with intensity measure nfi. Similarly X' — . . . ,X'j^_^} and 
y = {lY, . . . , Y^^} are independent Poisson point processes with intensity measure n/z'. 
The regularity property ensures that 

i({Xi,...,XjvJ,{ri,...,y^j) 
< L{{x[, . . . , x^j, {y/, . . . , r^}) + 2CAP Yl 1^.^^- + E ^y^^yj 



Taking expectations yields 



Afi N2 

EL{nn) < Ei( V) + 2CAPE I ^ ¥{X, ^ X'^) + ^ ¥{Yj ^ Y-) 



d\2. 



= EL{nfi')+ACAPmr{{{x,y)e{R'^f] x^y}). 

Optimizing the later term on the coupling tt yields the claimed inequality involving the total variation 
distance. □ 

Corollary 9. Assume that the functional L satisfies the regularity property (TZp). Let m > 0, Q C R'^ 
be measurable with positive Lebesgue measure and let f be a nonnegative locally integrable function on 
M.'^ . Let ct ~ Jq //vo1(Q) be the average value of f on Q. It holds 

EL{mflQ)<EL{malQ) + 2Cmdiam{Q)P / \f{x)~a\dx. 

Proof. We simply apply the total variation bound of the previous lemma with n = m f — mavol{Q), 
dii{x) = f{x)lQ{x)dx/ f and d^'{x) = lQ{x)dx/vo\{Q). Note that 

Jo \fi^) - "I '^^ 

dx — 



J in! VOI(Q) j / 



□ 
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2.1.3 Average is enough 

It is known since the works of Rhee and Talagrand that concentration inequalities often aUow to deduce 
ahnost sure convergence from convergence in average. This is the case in our general setting. 

Proposition 10. Let L be a bipartite functional on multisets of W^, satisfying the regularity property 
{Ti-p). Assume d > 2p > 0. Let ^ be a probability measure /i on M.'^ with J \x\"d^{x) < +oo. Consider 
independent variables {Xi)i>i and (Yi)i>i with distribution ^. 

If a > 2dp/{d— 2p) then the following convergence holds in probability: 



Moreover if a > 4:dp/{d—2p), the convergence happens almost surely, and if fi has bounded support, then 
it also holds in L'^ for any q>l. 

Proof. This is a simple consequence of Azuma's concentration inequality. It is convenient to Z{n) = 
(Ai, . . . , A„, Yi, . . . , Yn). Assume first that the support of fi is bounded and let A denote its diameter. 
By the regularity property, modifying one point changes the value of the functional by at most a constant: 

\L{Zi, . . . , Z2n) — L{Zi, . . . , Zi^i, Z'l, Zi+i, . . . , Z2n)\ < 2CA^. 

By conditional integration, we deduce that the following martingale difference: 

d, :- E(L(Z(n)) | Zi, . . . , Z,) - E(L(Z(n)) \Zi,..., Z,_i) 
is also bounded \di\ < 2CAP almost surely. Recall that Azuma's inequality states that 



>t\ <2e ^Eiiidiii 



Therefore, we obtain that 

|L({A,}r=i,{mi)-Ei({^Jr=i,{mi)| >i) <2e- 
and there is a number C (depending on A only) such that 



>t \ <2e 



-C't^n 



1-^ 



(2) 



When d > 2p, we may conclude by the Borel-Cantelli lemma. 

If fi is not assumed to be of bounded support, let 5 := max{|Zi|; i < 2n}. A conditioning argument 
allows to use the above method. Let s > and B{s) = {x; \x\ < s}. Given {S < s}, the variables 
{Ai, • • • , A„} and {Yi, ■ ■ ■ ,Yn} are mutually independent sequences with distribution ii\b{s)/ l^(.B{s)). 
Hence, applying ([2]) for p,\B(s)/ l^{B{s)) instead of ^ and 2s instead of A, for any i > 0, 



> t 



S <s \ < 2 exp 



Hence for (5 > to be chosen later, 

_ _ ^ / i({A,}r^i, {y.if^O EL({A,}r^i, {y.ir^O 



> t 



< V{S > n^) +2exp 



n^^ — ^s-t^ 



-tt' 



CpS^P 
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Since ¥{S > u) = I - (l - /i(B(s))2" < 2nfi{B{s)) < 2n{J we get that for some constant 

c and any 5 > 0, 



Un < cn^ ° + 2 exp 



Since a > 2dp/{d— 2p) we may choose 5 £ [2dp/{d— 2p),a], which ensures that the latter quantities 
tend to zero as n increases. This shows the convergence in probabihty to of 

1_£ 1_E 

n d n d 

If a > Adp/ {d — 2p) we may choose we may choose 6 G [2dp/{d — 2p), a/2] , which ensures that J2n ""'i 
+00. The Borel-CantcUi lemma yields the almost sure convergence to 0. □ 

2.2 Consequences of subadditivity 

We start with a very general statement, which is however not very precise when the measures do not 
have disjoint supports. 

Proposition 11. Let L satisfy (Sp). Let pi,fi2 be finite measures on R"^ with supports included in a set 
Q. Then 

EL{pi + ^2) < ELipi) + EL(/i2) + 2Cdiam(Q)P (^1 + ^/fii{Q) + y^jMOf) ■ 

Proof. Consider four independent Poisson point processes Xi,yi, X2,y2 such that for i e {1,2}, the 
intensity of Xi and of y.i is /i;. It is classical [5] that the random multiset Xi U X2 is a Poisson point 
process with intensity ^1 + /^2- Also, U 3^2 is an independent copy of the latter process. Applying the 
subadditivity property, 

LiXiuX2,yiuy2) 

< L{Xi,yi) + L{X2,y2) + Cdiam(Q)P (1 + |card(A'i) - card(3;i)| + 1 + |card(A'2) - card(3^2)|) ■ 
Since card(A'i) and card(3^i) are independent with Poisson law of parameter ^i{Q) (the total mass of /li), 

E|card(A;) - card(3^,)| < (^E(card(A'i) - card(yi))^) ' = V2var(card(A;)) === y/2fi,{Q). 

Hence, taking expectations in the former estimate leads to the claimed inequality. □ 

Partition techniques are essential in the probabilistic theory of Euclidean functionals. The next 
statement allows to apply them to bipartite functionals. In what follows, given a multiset X and a set 
P, we set X{P) := card(A' HP). If is a measure and / a nonnegative function, we write / • /i for the 
measure having density / with respect to fi. 

Proposition 12. Assume that the functional L satisfies (Sp). Consider a finite partition Q = Up^-pP 
of a subset o/R'^ and let v he a measure on with v{Q) < +00. Then 

EL{1q ■iy)<Y^ EL(lp • 1^) + 3Cdiam(Q)*' ^ V'^iP)- 

Pev pev 

Proof. Consider X,y two independent Poisson point processes with intensity i/. Note that A" n P is a 
Poisson point process with intensity Ip ■ ly, hence X{P) is a Poisson variable with parameter i^iP). We 
could apply the subadditivity property to {X n P)p^-p, {y n P)p^-p, which yields 

L{xr\Q,yr\Q)< ^ l{x r\P,yr\P) + Cdiam(Q)p + \x{p) - y{P)\). 

Pev pev 

Nevertheless, doing this gives a contribution at least Cdiam((5)^ to cells which do not intersect the 
multisets X,y. To avoid this rough estimate, we consider the cells which meet at least one of the 
multisets: 

■P := {P e V; X{P) + y{P) ^ 0}. 
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We get that 

L{xnQ,ynQ) < l{x n p,y n p) + Cdiam{Q)p"^ {i + \x{P) -y{p)\) 

Pev pGV 

< ^ LIA- n P,3; n P) + Cdiam(Of ^ lxiP)+y(P)M^ + \x{P) - y{P)\) 
Pev pev 

< L{xnP,ynP) + Cdiam{Q)pY.{ixiP)+yiP)^a + \x{P)-y{P)\). 

Pev pev 
Since X{P) and y{P) are independent Poisson variables with parameter z^(P), 



P{X{P)+y{P) 7^0) = l-e~2i^(^) andE\X{P)-y{P)\ < ^j2v(P). 
Hence, taking expectation and using the bound 1 — e^* < min(l,t) < \/<, 

EL(1q ■v)<Y EL(lp • v) + 2\/2Cdiani(Q)P ^ V'j^P). 

PeP pev 



□ 



The next statement deals with iterated partitions, which are very useful in the study of combinatorial 
optimisation problems, see e.g. [TS1[T1]. If 7-" is a partition, we set diam(7') = maxpg-p diam(P) (the 
maximal diameter of its cells). 

Corollary 13. Assume that the functional L satisfies (Sp). Let Q C M'* and Qi, . . . , Qk he a sequence 
of finer and finer finite partitions of Q. Let v he a measure on with i'{Q) < +oo. Then 

k 

EL(1q ■v)<Y. ^^(1? ■^) + ^CY. diam(Q,_i)P ^ ^^9), 

where hy convention Qo = {Q} is the trivial partition. 

Proof. We start with applying Proposition [T^ to the partition Qi of Q: 

EL(1q ■v)<Y^ EP(1, • v) + 3Cdiam(Qo)^ ^ 

qeQi qeQi 

Next for each g G Qi we apply the proposition again for the partition of q induced by Q2 and iterate 
the process k — 2 times. □ 

3 Uniform cube samples 

We introduce a specific notation for n G (0, +00), 

L{n) := EL(?il[o.i]ti) . 

We point out the following easy consequence of the homogeneity properties of Poisson point processes. 

Lemma 14. If L satisfies the homogeneity property (Hp) then for all a G M'', p > Q and n > 

EL(nl,+[o,p].) ^pPL{np'^). 

The following theorem is obtained by adapting to our abstract setting the line of reasoning of Boutet 
de Monvcl and Martin in the paper [3] which was devoted to the bipartite matching: 

Theorem 15. Let d > 2p be an integer. Let L be a bipartite functional on R'^ satisfying the properties 
{T-Lp), (Ji-p) and (Sp). Then there exists /3l > such that 

lim4^=/3.. 
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Proof. Let m > 1 be an integer. Let if G N such that 2^ < m < 2^+^ Set Qo = [0,0]^^ where 
a :— 2^'+^/m > 1. Let Qq = {Qo}- We consider a sequence of finer and finer partitions Qj, j > 1 where 
Qj is a partition of Qq into 2^'^ cubes of size a2~^ (throughout the paper, this means that the interior 
of the ceUs are open cubes of such size, while their closure is a closed cube of the same size. We do not 
describe precisely how the points in the boundaries of the cubes are partitioned, since it is not relevent 
for the argument). One often says that Qj, j > 1 is a sequence of dyadic partitions of Qq. 

A direct application of Corollary 1131 for the partitions Qi, . . . , Qk+i and the measure nlij^ ^d(x) dx 
gives 



Z(n) = EX(nl[o,i].)< ^ EL(nl,n[o,i]'') + 3C ^ diam(Q,_i)P ^ Vol(g n [0, 1]^). 

Note that Qk+i is a partition into cubes of size 1/m, so that its intersection with [0,1]'' induces an 
(essential) partition of the unit cube into m'^ cubes of side-length 1/m. Hence, in the first sum, there 
are m'^ terms which are equal, thanks to translation invariance and Lemma 1141 to Eij(nl[Q ,„-i]<i) = 
m~PL{nm~'^). The remaining terms of the first sum vanish. In order to deal with the second sum of the 
above estimate, we simply use the fact that Qj contains 2^'^ cubical cells of size a2^^ = 2^+^^-' /m < 2^~^ . 
Hence their indidual volumes are at most 2'^^^^^\ These observations allow to rewrite the above estimate 
as 

K+l 

L{n) < m'^-PL{nm-'^) + 3C ^ diam([0, 2'^-^]y2^'^\/n2di^-^) 

A'+l 

= m'^-PL{nm-'^) + 3Cv^diam([0, 1]'')^ ^ 2f(2-J)+lO"+i). 

i=i 

Hence, there is a number D depending only on p, d and C such that 

L{n) < ni'^-PL{nm-'^) + D^/^2'^'^i'P^ < m'^-PL[nm-'^) + Dy/^m^-P. 
Let t > 0. Setting, n — m'^t'^ and f{u) = L{u'^) / u'^~p , the latter inequality reads as 

fimt)<fit) + DtP-i, 

and is valid for all f > and m G N*. Since / is continuous (Proposition [6] shows that u H> L{u) is 
Lipschitz) and limt^+oo tP^^ = 0, it follows that limt^^+oo f{t) exists (we refer to [3] for details). □ 

Remark 4. The above constant /3l is positive as soon as L satisfies the following natural condition: for 
ah xi,...,a;„,yi,...y„ in R'', L{{xi, . . . ,x„},{yi, . . . ,yn}) > cJ2i<iist(xi, {yi, . . . ,yn})P ■ To see this, 
one combines Proposition [5] and the lower estimate given in |16) . 

4 Upper bounds, upper limits 

4.1 A general upper bound 

Lemma 16. Let d > 2p and let L be a bipartite functional satisfying {Sp), (TZp) and L(0,0) = 0. Then 
there exists a constant D such that, for all finite measures v, 

EL{v) < D diam{Q)P min {iy{Q),v{Qy-^) , 

where Q contains the support of v. 

Proof. Thanks to corollary[71 it is enough to deal with the case v{Q) > 2"^ (or any other positive number). 
First note that we may assume that Q is a cube (given a set of diameter A, one can find a cube containing 
it, with diameter no more than c times A where c only depends on the norm). We consider a sequence 
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of dyadic partitions of Q, {'Pe)e>o, where for £ £N,Ve divides Q into 2^'^ cubes of side-length 2 ^ times 
the one of Q. Let S N* to be chosen later. By Corollary [T31 we have the following estimate 

k 

EL(z^) < ^ EL(lp ■iy)+3Cj2 (2-^+Miam(Q))'' ^ ^MP). (3) 

P£Vk 1=1 PeVe 

Thanks to Corollary [71 the first term of the right-hand side of ([3]) is at most 

^ C(2-'^diam(g))^i/(P) =C2-*^f(diam(Q))''zy(Q). 

Pen 

By the Cauchy-Schwarz inequality 

Hence the second term of the right-hand side of ([3]) is at most 

3C(2diam(Q))^ J2 2'^^'") ^J^Aff) < C'2^(^-p) {diam{Q)Y . 
1=1 

This leads to 

¥.L{v) < {dmmiQ)y (c2-''Piy{Q) + C'2''ii~p) 
Choosing k — \_2 log2 '^iQ)\ — 1 completes the proof. □ 

4.2 The upper limit for densities 

Theorem 17. Let d > 2p. Let L he a bipartite functional on satisfying the properties {Hp), O^p), 
(Sp). Let f : — )■ IR+ be an integrahle function with bounded support. Then 

hm sup — < Pl f 
where jS^ is the constant appearing in Theorem \15\ 

Proof. By a scaling argument, we may assume that the support of / is included in [0, 1]'' and J f = 1 
(the case J / = is trivial). We consider a sequence of dyadic partitions (7-'f)fgN of [0, 1]'': for £ £N,Ve 
divides [0, 1]'' into 2^'^ cubes of side-length 2~^. Let A; G N* to be chosen later. Corollary [TS] gives 



Ei(n/)< J2 mnflp) + SCY,{2-'+'dmm{[0,lf)y ^ Jn f f. (4) 



PeVk 1=1 PeVe 

By the Cauchy-Schwarz inequality 



Hence the second term of the right-hand side of is at most 

3C(2diam([0,l]^))^V^^2<^-?') < c^nh^i-p) . 

1=1 

Let ap be the average of / on P, then applying Corollary [3] to the first terms of ^ leads to 

EL{nf)< (EL{naplp) + 2CndiSim{P)P f \f - ap\\ + Cdnh'^^^'P) . 

PeVk ^ J P J 
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Each P in the sum is a square of side length 2 hence using homogeneity (see Lemma [H)) 

Ei(n/)< ^ (2~^PM{nap2-^'^) +nc'^2-''P f \f - ap\) + Cdnh''(i~p) . (5) 
PePfc V Jp / 

Let us recast this incquahty with more convenient notation. We set g{t) = L{t)/t^~P^'^ and we define 
the piecewise constant function 

, ^ , ^ r„ f(x) dx 

P£Vk PeVk ^ ' 

It is plain that J fk = J f < +00. Moreover, by Lcbcsguc's theorem, limfe^oo fk = f holds for almost 
every point x. Inequality (jSj amounts to 

n Jp J 



PeVk 

g(n2-'='^/fe)/r^ + c',n^- 2-"^ j \f - M + Crf(n^2-'^y-*. 

If there exists fcg such that / — fka then we easily get the claim by setting k = ko and letting n go to 
infinity (since g is bounded and converges to (5^ at infinity, see Lemma 1161 and Theorem II 5p. On the 
other hand, if fk never coincides almost surely with /, we use a sequence of numbers k{n) g N such that 

limA:(n) = +00, limn^2-'=(") = +00 and limn32-'=(") ( f \f ~ fk(n)\] " = 0. (6) 

n n n \J J 

Assuming its existence, the claim follows easily: applying the inequality for k = kin) and taking upper 
limits gives 

limsup?^^ < hmsup / 5(r^2-'=(")'^/,(„))4\;f . 

Since Mmf^n) = f a.e., it is easy to see that the limit of the latter integral is /3l / /^~^: first the 
integrand converges almost everywhere to (if fix) = this follows from the boundedness of g; if 

f{x) 7^ then the argument of g is going to infinity). Secondly, the sequence of integrands is supported 
on the unit cube and is uniformly integrable since 

(,9(n2-^(")''/fc(„))/r^)'^ < (sup5)^ J fkin) = (supg)^ J f < 

It remains to establish the existence of a sequence of integers (fc(n))„ satisfying Note that since 
fk ^ 0, J fk = / / = 1 E^nd a.e. lim/fc = /, it follows from Scheffe's lemma that limfc J \ f — fk\ = 0. 
Hence (p{k) = (sup^^j, / |/ — fjl)^"^^^ is non-decreasing with an infinite limit. We derive the existence of 
a sequence with the following stronger properties 

n, Ti 
limfc(n)=+oo, lini^^^ = +oo and 1™ = » (7) 



as follows. Set 7 = 2'*. Since 7*^^/ ~ 1) is increasing with infinite limit 



[7 VRo), +00) = Ufc>i [7^^ V^(fc-i), 7'+' V^) ■ 



For n > 7^ (p{0), we define k{n) as the integer such that 



7'=(") ^/^(Mn)^ <n< 7*^(")+i V^{k{n)). 

This defines a non-decreasing sequence. It is clear from the above strict inequality that lim„ k(n) = +00. 
Hence n7~'^^"-' > ip{k{n) — 1) tends to infinity at infinity. Eventually n/{-/'^'^"^ip{k{n)) < ^ / ip{k{n)) 
tends to zero as required. The proof is therefore complete. □ 
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4.3 Purely singular measures 

Lemma 18. Let d > 2p. Let L be a bipartite functional on R'^ with properties (TZp) and {Sp). Let fi he 
a finite singular measure on having a bounded support. Then 

lim — = 0. 

Proof. Let Q be a cube which contains the support of ^. We consider a sequence of dyadic partitions 
of Q, {'Pe)e£n. For £ € N, Ve divides Q into 2^'^ cubes of side length times the one of Q. As in the 
proof of Lemma [T51 a direct appUcation of Corollary [T^ gives for € N*: 

fe 

EL(n^) < Ei(nlp • /i) + 3C^ (2-''+Miam(g))^ ^ V^niP). (8) 

PeVk 1=1 PeVi 

The terms of the first sum arc estimated again thanks to the easy bound of Corollary [T] since each P in 
Vk is a cube of side length 2~'' times the one of Q, it holds 

EL{nlp-n)< Y C{2-''diam{Q)Ynfi{P) ^ Cp^Q2-''Pn\fi\. 
Here |^| is the total mass of ^. We rewrite the second term in ([8]) in terms of the function 

PeVt ^ ' 

where A stands for Lebcsgue's measure. Since A(P) = 2~'^^\(ff). we get that 

k 



A(P) 



EL(n/i) < Cp,Q2-'^■Pn|/i|+3C(2diam(Q))''^^^2-^P ^ 2^A(g)-U(P)i 

1=1 P^Vt. 
k 

= Cp,Q2-'^-Pn|M|+3C(2diam(Q))^A(0)-^V^^2^(*-'') j^,. 

1=1 •' 

By the differentiability theorem, for Lebcsgue-almost every x, ge{x) tends to zero when i tends to infinity 
(since fj, is singular with respect to Lcbesgue's measure). Moreover, gi is supported on the unit cube and 
Jiy/gl)^ = J gt = < +0O. Hence the sequence of functions y/gj is uniformly integrable and we can 
conclude that \imi=foo J ^/gl = 0. By Cesaro's theorem, the sequence 

also converges to zero, using here that d > 2p. By an obvious upper bound of the latter denominator, 
we obtain that there exists a number c which does not depend on {k,n) (but depends on C\p,d,Q, 
such that for all A; > 1 

EL{nn) < c (n2-''P + y^i2''i^~p) Sk 



where Sk 1^ and limfc Ek = 0. We may also assume that (sk) is non-increasing (the inequality remains 
valid if one replaces by supj>j. Sj). It remains to choose k in terms of n in a proper way. Define 

-1 

Obviously lim„ ip{n) ~ +oo. For n large enough, define k(n) > 1 as the unique integer such that 

2'^(") < n^ifin) < 2'=(")+i. 
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Setting k = k{n), our estimate on the cost of the optimal matching yields 



It is easy to check that the right hand side tends to zero as n tends to infinity. Indeed, lim„ ip{n) = +00, 
hence for n large enough 

k{n) > [log2 (?i^(^(n)/2)J ^ ^ 



> 



Since the sequence (sk) is non-increasing, it follows that 

tends to zero when n — > 00. The proof is therefore complete. □ 
4.4 General upper limits 

The first statement of Theorem [3] is a consequence of Propositions [51 [TUl and the following result. 

Theorem 19. Let d > 2p > Q. Let L be a bipartite functional on M/^ with the properties {Hp), {Ttp) and 
(Sp). Consider a finite measure pL on such that there exists a > -jz^ with 



\x\°'dpi{x) < +00. 

Let f be a density function for the absolutely continuous part of fi, then 

EL(nu) ^ f p 
limsup ;_r <fa / /'~^- (9) 

Remark 5. Observe that the hypotheses ensure the finiteness of J f^^^ ■ Indeed Holder's inequality gives 

f-^<(f (i + \xnf{x)dx) Uf {i + \xn'-i 

\Jr'' J \Js.d 

where the latter integral converges since a > > -jiz^- 

Proof. Assume first that fi has a bounded support. Write yU = /.iac + Ms where is the singular part 
and dfj.ac{x) = f{x) dx. Applying Proposition 1 1 II to ^ac and pis, dividing by v}~'pI'^ , passing to the limit 
and using Theorem [T7| and Lemma [18] gives 

limsup — _p < Inn sup p + hmsup ' p < h f '' ■ 



Hence the theorem is established for measures with bounded supports. 

Now, let us consider the general case. Let B{t) = {x e M'' : < t}. Let Aq = B{2) and for integer 
£ > 1, Ae = B(2^+i)\B(20. Now, let X = {Xi,--- ,XnA, 3^ = {>!,••• , >V J be two independent 
Poisson process of intensity ufx, and T = max{|Z| : Z e X uy}. Applying the subadditivity property 
like in the proof of Proposition [T^ we obtain 

L{X,y) < Y.LiXnA,,ynAe) + CTPY,^xiA,)+yiA,)M^ + \X{A,)-yiAi)\). (10) 

Note that the above sums have only finitely many non-zero terms, since fj, is finite. We first deal with 
the first sum in the above inequality. By Fubini's Theorem, 

L{xnAi,ynAe) _y^^ L{xnAi,ynAe) 
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Applying ^ to the compactly supported measure fj.^^^ for every integer £ gives 

hmsupE^ — p -<Pl f 

n d J Ai 

By Lemma [TOl for some constant Cd, 

,L{X r^ At,y r\ At) 



(11) 



E- 



<cpJp^l{A,f-^d, 



From Markov inequality, with 7ti„ = J \x\'^d^{x), 

Thus, since a > 2pd/{d — 2p) > dp/{d — p), the series J2e '2'''^fJ-iAi)^^^ is convergent. We may then apply 
the dominated convergence theorem, we get from that 

hmsupE) — £ </3l// " ■ 

^ 3 J 

For the expectation of the second term on the right hand side of pIT]) . we use Cauchy- Schwartz inequality, 



E 



{l + \X{Ae)-y{A,)\) 



< J2 VW^WH^xiA,)+yiA,)M^ + \XiAe) - yiA,)\Y) 

< v^VfifT^^ ^nx{Ae) + y{Ae) + 0) + E[|A'(A,) - y{A,)Y\ 



£>0 



< 



£>0 



where we have used 1 — e " < u. As above, Markov inequality leads to 



Eventually we apply Lemma 2] with 7 := 2p < 2pd/ (d — 2) < a to upper bound E[T^^']. We get that for 
some constant c > and all n > 0, 



n-i+SE 



TPY.lxiA,)+yiA,)M^ + \XiA,) - yiA,)\) 



< cn 2 + d " 



Since a > 2dp/{d— 2p), the later and former terms tend to zero as n tends to infinity. The upper bound 
© is proved. □ 

5 Examples of bipartite functionals 

The minimal bipartite matching is an instance of a bipartite Euclidean functional Mi{X,y) over the 
multisets X = {Xi, . . . , Xn} and y = {Fi,...,F„}. We may mention at least two other interesting 
examples: the bipartite traveling salesperson problem over X and y is the shortest cycle on the multiset 
X U y such that the image of X is y. Similarly, the bipartite minimal spanning tree is the minimal 
edge-length spanning tree on A" U 3^ with no edge between two elements of X or two elements of y. 
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5.1 Minimal bipartite matching 

Fix p > 0. Given two multi-subsets of Mf^ with the same cardinahty, X = {Xi, . . . , Xn} and y = 
{Yi, . . . ,Yn}, the p-cost of the minimal bipartite matching of X and y is defined as 

n 

Mp{X,y)^ min V|X,-y,(,)|P, 

i=l 

where the minimum runs over all permutations of {1, . . . , ri}. It is useful to extend the definition to 
sets of different cardinalities, by matching as many points as possible: if X ~ {Xi, . . . ^Xm} and y = 
{Yi, . . . , Yn} and m < n then 



MpiX,y)^mmJ2\X^-Ya 



where the minimum runs over all injective maps from {1,...,™} to When n < m the 

symmetric definition is chosen Mp{X,y) := Mp{y,X). 

The bipartite functional Mp is obviously homogeneous of degree p, i.e. it satisfies (Tip). The next 
lemma asserts that it is also verifies the subadditivity property {Sp). In the case p = 1, this is the 
starting point of the paper 



Lemma 20. For any p > 0, the functional Mp satisfies property {Sp) with constant C = 1/2. More 
precisely, if Xi, . . . , Xk and 3^i, . . . ,3^fc are multisets in a bounded subset Q C M.'^, then 

Mp( y X,, \Jy^)<Y. AMX^^y^) + 2 ^ IcardCA-,) - card(3;,)l- 

i—l i—1 i—1 i—1 

Proof. It is enough to upper estimate of the cost of a particular matching of IJ^^^ Xi and IJ^^^ ^i- ^® 
build a matching of these multisets as follows. For each i we choose the optimal matching of Xi and 
yi. The overall cost is J^i Mp{Xi,yi), but we have left J^i |card(A:i) — card(3^i)| points unmatched (the 
number of excess points). Among these points, the less numerous species (there are two species: points 
from Xi's, and points from J^^'s) has cardinality at most 5 |card(A'i) — card(3^i)|. To complete the 
definition of the matching, we have to match all the points of this species in the minority. We do this in 
an arbitrary manner and simply upper bound the distance between matched points by the diameter of 
Q. □ 

The regularity property is established next. 



Lemma 21. For any p > 0, the functional Mp satisfies property [TZp I with constant C = 1. 

Proof. Let X , Xi, X2,y,yi,y2 be finite multisets contained in Q = 5(1/2). Denote by x,a;i,a;2,y,yi,y2 
the cardinalities of the multisets and a A 6 for min(a, b). We start with an optimal matching for Mp{X D 
'^2,y n 3^2)- It comprises {x + X2) /\ {y + 2/2) edges. We remove the ones which have a vertex in X2 or 
in 3^2- There are at most X2 + 2/2 of them, so we are left with at least ((a; + X2) A (y + 2/2) ~ X2 ~ 2/2)^ 
edges connecting points of X to points of y. We want to use this partial matching in order to build a 
(suboptimal) matching of X C] Xi and 3^ n J^i. This requires to have globally {x + xi) A {y + yi) edges. 
Hence we need to add at most 

{x + xi) A{y + yi) - {{x + X2) A (y + ^2) - X2 - ^2)^ 

new edges. We do this in an arbitrary way, and simply upper bound their length by the diameter of Q. 
To prove the claim it is therefore sufficient to prove the following inequalities for non-negative numbers: 

{x + xi) A{y + yi) - {{x + X2) A {y + 1/2) - 2^2 - y2)_^_ < xi + X2 + yi + ?/2- (12) 

This is obviously equivalent to 

x + xi < xi + X2 + yi + y2 + {{x + X2) A{y + 2/2) - X2 - ^2)^ 
or 2/ + 2/1 < 2^1 + 2^2 + 2/1 + 2/2 + {{x + X2) A (2/ + 2/2) - 2:2 - 2/2) 
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After simplification, and noting that yi > appears only on the right-hand side of the first inequation 
(and the same for xi in the second one), it is enough to show that 

X Ay < X2 +y2 + {{x + X2) A {y + ^2) - X2 - j/2)_^- 

This is obvious, as by definition of the positive part, x A y < X2 + y2 + {{x A y) — X2 — 2/2) + - D 

5.2 Bipartite traveling salesperson tour 

Fix p > 0. Given two multi-subsets of M"^ with the same cardinality, X ~ {Xi, . . . , Xn} and y = 
{Yi, . . . , Yn}, the p-cost of the minimal bipartite traveling salesperson tour of {X ,y) is defined as 

Tp{X,y) = min I^-W " ^-'Wl'' + E " ^-C^+dT + " ^-(dT, 

where the minimum runs over all pairs of permutations of {1, . . . , n}. We extend the definition to sets 
of different cardinalities, by completing the longest possible bipartite tour : ii X = {Xi, . . . ,Xm} and 
y = {Yi , . . . ,Yn} and m <n then 

m m — 1 

Tp{x,y) = min V - + V - + I^^.'m - ^.(dT 

where the minimum runs over all pairs (ct, a'), with a € Sm and a' is an injective maps from {1, . . . , m} 
to {l,...,n}. When n <m the symmetric definition is chosen Tp{X ^y) := Tp{y,X). This traveling 
salesperson functional is an instance of a larger class of functionals that wc now describe. 

5.3 Euclidean combinatorial optimization over bipartite graphs 

For integers ?ti, 71, we define [n] = {1, ■ ■ ■ n\ and [n\m = {m-\- 1, • ■ • , m + n\. Let S„ be the set of bipartite 
graphs with common vertex set ([n], : if G G S„, the edge set of G is contained is the set of pairs 
{i,7i-t- j}, with i,j e [n]. 

We should introduce some graph definitions. If Gi e Bn and G2 € ;Bm we define Gi -I- G2 as the 
graph in Bn+m obtained by the following rule : if {i, n + j} is an edge of Gi then {z, n + m + j} is an edge 
of Gi + G2, and if {i, m + j} is an edge of G2 then {?? + i, 2n + m + j} is an edge of Gi 4- G2. Finally, if 
G e Bn+m; the restriction G' of G to Bn is the element of Bn defined by the following construction rule: 
if {i, ji + m + j} is an edge of G and (i, j) € [n]^ then add {?, n -|- j} as an edge of G'. 

We consider a collection of subsets Gn C /S„ with the following properties, there exist constants 
i^Ot 1^ 1^ 1 such that for all integers n, m, 

(Al) (not empty) If 71 > kq, Gn is not empty. 

(A2) (isomorphism) If G G Gn and G' e S„ is isomorphic to G then G' E Gn- 
(A3) (bounded degree) If G G the degree of any vertex is at most k. 

(A4) (merging) li G G Gn and G' G f/m, there exists G" G Gn+m such that G -I- G' and G" have all but 
at most K edges in common. For 1 < m < kqj it also holds if G' is the empty graph of Bm. 

(A5) (restriction) Let G G Gn and kq + \ <n and G' be the restriction of G to Bn-i- Then there exists 
G" G Gn-i such that G' and G" have all but at most n edges in common. 

If I A- 1 = \y\ = n, we define 

L(A',3^)=min ^ 

" {i.'j)e[n]'^:{i,n+j}eG 

With the convention that the minimum over an empty set is 0. Note that the isomorphism property 
implies that L[X,y) = L{y,X). If m = \X\ < \y\ = n, we define 

L(A-,3^)= min ^ \X, - Y^^.^l" , (13) 

{i,])el7n]^:{i,7n+j}eG 
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where the miminum runs over all pairs (G, a), G G Qm and a is an injective maps from {1, . . . , m} to 
{!,..., n}. When n < m the symmetric definition is chosen L{X,y) := L{y, X). 

The case of bipartite matchings is recovered by choosing Qn as the set of graphs in i3„ where all 
vertices have degree 1. We then have kq ~ \ and Qn satisfies the merging property with k = 0. It also 
satisfies the restriction property with k = 1. The case of the traveling salesperson tour is obtained by 
choosing C/„ as the set of connected graphs in Bn where all vertices have degree 2, this set is non-empty 
for n > kq = 2. Also this set Qn satisfies the merging property with k = 4 (as can be checked by edge 
switching). The restriction property follows by merging strings into a cycle. 

For the minimal bipartite spanning tree, we choose Qn as the set of connected trees of [2n] in Bn- 
It satisfies the restriction property and the merging property with k — 1. For this choice, however, the 
maximal degree is not bounded uniformly in n. We could impose artificially this condition by defining 
Qn as the set of connected graphs in Bn with maximal degree bounded by k > 2. We would then get the 
minimal bipartite spanning tree with maximal degree bounded by k. It is not hard to verify that the 
corresponding functional satisfies all the above properties. 

Another interesting example is the following. Fix an integer r > 2. Recall that a graph is r-regular 
if the degree of all its vertices is equal to r. We may define Qn as the set of r-regular connected graphs 
in Bn- This set is not empty for n > kq ^ r- It satisfies the first part of the merging property (A4) 
with K = 4. Indeed, consider two r-regular graphs G, G', and take any edge e = {x,y} G G and 
e' = {a;', y'} £ G'. The merging property holds with G", the graph obtained from G 4- G' by switching 
(e, e') in {{x,y'},{x',y})- Up to increasing the value of k, the second part of the merging property is 
also satisfied. Indeed, if n is large enough, it is possible to find rm < vkq = edges ei, • • • , e^m in 
G with no-adjacent vertices. Now, in G", wc add m points from each species, and replace the edge 
Gri+q = {x, n + J/}, 1 < i < 777 , < g < r, by two edges : one between x and the i-th point of the second 
species, and one between y and the i-th point of the first species. G" is then a connected r-regular graph 
in Bn+m with all but at most 2r^ edges in common with G. Hence, by taking k large enough, the second 
part of the merging property holds. 

Checking the restriction property (A5) for r-regular graphs requires a little more care. Let r = 
kq + 1 < n and consider the restriction Gi of G G Bn to Bn-i- Our goal is to show that by modifying 
a small number of edges of Gi, one can obtained a connected r-regular bipartite graph on Bn-i- We 
first explain how to turn Gi into a possibly non-connected r-regular graph. Let us observe that Gi was 
obtained from G by deleting one vertex of each spieces and the edges to which these points belong. Hence 
Gi has vertices of degree r, and vertices of degree r — 1 (r blue and r red vertices if the removed points 
did not share an edge, only r — 1 points of each spieces if the removed points shared an edge). In any 
case Gi has at most 2r connected components and r vertives of each color with one edge missing. The 
simplest way to turn Gi into a r regular graph is to connect each blue vertex missing an edge with a red 
vertex missing an edge. However this is not always possible as these vertices may already be neighbours 
in Gi and we do not allow multiple edges. However given a red vertex vr and a blue vertex vb of degree 
r — 1 and provided n — \ > 2r^ there exists a vertex v in Gi which is at graph distance at least 3 from 
Vb and vr. Then open up an edge to which v belongs and connect its end-points to vr and vb while 
respecting the bipartite structure. In the new graph vb and vr have degree r. Repeating this operation 
no more than r times turns Gi into a r regular graphs with at most as many connected components 
(and the initial and the final graph differ by at most 3r edges). Next we apply the merge operation 
at most 2r — 1 times in order to glue together the connected componented (this leads to modifying at 
most 4(2r — 1) edges. As a conclusion, provided we choose kq > 2r^, the restriction property holds for 
K = llr. 

We now come back to the general case. From the definition, it is clear that L satisfies the property 
\H.p\- We are going to check that it also satisfies properties and (|7^pp. 



Lemma 22. Assume (A1-A4-)- For any p > 0, the functional L satisfies property [Sp) with constant 
C= (3 + ko)k/2. 

Proof- The proof of is an extension of the proof of Lemma 1201 We can assume without loss of generality 
k > 2. Let Xi, . . . ,Xk and J^i, . . . ,3^fe be multisets in Q = _B(l/2). For case of notation, let Xi — \Xi\, 
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Vi — \yi\ and n = Y^^^i A ^i^i Vi- If n < /^o, then from the bounded degree property (A3), 



i=l 1=1 



If n > kq, it is enough to upper bound the cost for an element G in C/„. For each 1 < i < k, ii 
rii = Xi A Ui > kq, we consider the element Gi in C/„. which reaches the minimum cost of L{Xi,yi). From 
the merging property (A4), there exists G' in Qj2 whose total cost is at most 

i 

It remains at most kq + \xi — yi\ vertices that have been left aside. The less numerous species has 
cardinal mo < m = (^,- kq + \xi — yi\)/2. If mo > kq, from the non-empty property (Al), there exists 
a graph G" G Gmo that minimizes the cost of the vertices that have been left aside. From the merging 
and bounded degree properties, we get 

k k 

L{\jx,,\jy,^<L' + n + Km<Y, W,y,) + I ^ (3 + + - y^\) ■ 

i—1 i—1 i i 

If mo < kq, wc apply to G' the merging property with the empty graph : there exists an element G in 
Qn whose total cost is at most 

k k 

l[ U X,, \Jy^) <L' + K<Y, L{x,,y,) + [k + 

i—1 i—1 i 

Wc have proved that property \Sp^ is satisfied for C = (3 + ko)k/2. □ 



Lemma 23. Assume (A1-A5). For any p > 0, the functional L satisfies property [TZp) with constant 
G ~ G{k, kq). 

Proof. Let X,Xi,X2,y,yi,y2 be finite multisets contained in _B(l/2) — Q. Denote by x,xi,X2,y,yi,y2 
the cardinalities of the multisets. As a first step, let us prove that 

L{X UXi,yu yi) < L{X,y) + + y^). (14) 

By induction, it is enough to deal with the cases (xi, j/i) = (1,0) and (xi, j/i) ~ (0, 1). Because of our 
symmetry assumption, our task is to prove that 

L{X\j{a},y)<L{X,y) + G. (15) 

If card(y) < card(A'), then the latter is obvious: choose an optimal graph for L{X ,y) and use it to 
upper estimate L{X U {a}, 3^). Assume on the contrary that card(3^) > card(A') + 1. Then there exists 
y cy with card(3^') = card(<Y) and L{X,y'). Let & G 3^ \ Y . In order to establish (HI]), it is enough 
to show that 

L{X U {a}, y U {6}) < L{X, y ) + C, 
but this is just an instance of the subadditivity property. Hence (jl4p is established. 
In order to prove the regularity property, it remains to show that 

L{X,y)<L{X\JX2,y\jy2) + G{x2+y2)- (16) 

Again, using induction and symmetry, it is sufficient to establish 

L{X,y)<L{X\j{a},y) + G. (17) 

If card(A') A card(cy) < ko, then by the bounded degree property L{X,y) < KKodia,m{Q)P and we are 
done. Assume next that card(A'), card(3^) > kq. Let us consider an optimal graph for L{X U {a},y). 
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If a is not a vertex of this graph (which forces card(A') > card(3^)) then one can use the same graph 
to upper estimate L{X,y) and obtain (jl7p . Assume on the contrary that a is a vertex of this optimal 
graph. Let us distinguish two cases: if card(A') > card(3^), then in the optimal graph for L{X U {a},y), 
at least a point b E X is not used. Consider the isomorphic graph obtained by replacing a by 6 while the 
other points remain fixed (this leads to the deformation of the edges out of a. There are at most k of 
them by the bounded degree assumption). This graph can be used to upper estimate L{X,y), and gives 

LiX.y) < L{XU{a},y) + Kdmm{Q)P. 

The second case is when a is used but card(<Y) + 1 < card(3^). Actually, the optimal graph for L{X U 
{a}, y) uses all the points of A'U {a} and of a subset of same cardinality y' C y. Choose an clement b in 
y. Theny = y'\{b} has the same cardinality as X. Obviously L(A' U {a}, 3^) = L{X U {a},y" U {b}). 
Consider the corresponding optimal bipartite graph. By the restriction property, if we erase a and b 
and their edges, we obtain a bipartite graph on {X,y") which differs from an admissible graph of our 
optimization problem by at most k edges. Using this new graphs yields 

LiX,y) < Kdmm{Q)P + L{X U {a},y" U {b}) ^ Kdiam{Qf + L{X U {a},y). 

This concludes the proof. □ 



6 Lower bounds, lower limits 
6.1 Uniform distribution on a set 

In order to motivate the sequel, we start with the simple case where / is an indicator function. The 
lower bound is then a direct consequence of Theorem [15] and Theorem 1171 

Theorem 24. Let d > 2p > 0. Let L be a bipartite functional on R'^ satisfying the properties (Hp), 
{Ti-p), {Sp). Let J7 C M'' be a bounded set with positive Lebesgue measure. Then 

hm 5^^ = /3,Vol(17). 

Proof. Theorem [T7l gives directly limsupEL(nlo)/?i^^^ < /3LVol(f2). By translation and dilation in- 
variance, we may assume without loss of generality that il C [0,1]**. Let ilc := [0,1]'' \ fi. Applying 
Proposition [T^ for the partition [0, 1]'^ = flU flc, gives after division by n^^P^'^ 



iL{nl,ni]d) EL(nln) ELlnln) ^ p i / i i 

' ^ ^ < ^3^+3Cdiam([0,l]'^)n§-5(Vol(r!)2 +Vol(r!,)2 



Since d > 2p, letting n go to infinity gives 

. EL(nlo) ^ EL(nl[o,i].) EL(nloJ 
hm mt — ^ — > Imi —3 lim sup ——^ 

where we have used Theorem [T51 for the limit and Theorem [T7] for the upper limit. □ 

The argument of the previous proof relies on the fact that the quantity limn^~^'/'^EL(nlsi) = 
I3lYo\{VL) is in a sense additive in i7. This line of reasoning does not pass to functions since / 1— >■ / f^^P/''- 
is additive only for functions with disjoint supports. The lower limit result requires more work for general 
densities. 



6.2 Lower limits for matchings 

In order to establish a tight estimate on the lower limit, it is natural to try and reverse the partition 
inequality given in Proposition 1121 This is usually more difficult and there does not exist a general 
method to perform this lower bound. We shall first restrict our attention to the case of the matching 
functional Mp with p > 0, we define in this subsection 

L = AL. 
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6.2.1 Boundary functional 

Given a matching on tlic unit cube, one needs to infer from it matchings on the subcubes of a dyadic 
partition and to control the corresponding costs. The main difficulty comes from the points of a subcube 
that are matched to points of another subcube. In other words some links of the optimal matching cross 
the boundaries of the cells. As in the book by Yukich [12], a modified notion of the cost of a matching is 
used in order to control the effects of the boundary of the cells of a partition. Our argument is however 
more involved, since the good bound ((T]) used by Yukich is not available for the bipartite matching. 

We define 

q = 2P-^ A 1. 

Let S C M.'^ and e > 0. Given multisets X — {Xi, . . . , m} and y = {Yi, . . . , y„} included in S we 
define the penalized boundary-matching cost as follows 

LdS,eiXl, . . . , Xm;Yi, . . . ,Yn) (18) 

' '"^ [isA ieA<= jeB'-- J 

where the minimum runs over all choices of subsets A C {l,...,m}, B C {l,...,ri} with the same 
cardinality and all bijective maps a : A ^ B. When e = we simply write Lqs- Notice that in our 
definition, and contrary to the definition of optimal matching, all points are matched even if m 7^ n. If 
X and y are independent Poisson point processes with intensity ly supported in S and with finite total 
mass, we write LQs,e{y) for the random variable L9s_e(A', 3^). 

The main interest of the notion of boundary matching is that it allows to bound from below the 
matching cost on a large set in terms of contributions on cells of a partition. The following Lemma 
establishes a superadditive property of Lgs a-nd it can be viewed as a counterpart to the upper bound 
provided by Proposition [T2l 

Lemma 25. Assume L = A/p. Let v he a finite measure on and consider a partition Q = UpgpP of 
a subset ofR^. Then 

diam(Q)Py^2i^(R'^) + EL{iy) > ELagilQ ■ v) > ELap{lp ■ v). 

Per 

Proof. Let X = {^1, . • . , X^}, y = {Yi, . . . , Yn} be multisets included in Q and X' ~ {Xm+i, ■ ■ ■ , Xm+m'}, 
y = {Fri+i, • ■ • , Yn+n'} bc multisets included in Q"^. By considering an optimal matching oi XUX' and 
3^ U 3^', we have the lower bound 

diam(Q)P|m + m' - n - n'\ + L{X U X' ,y U y') > LaQ{X,y). 

Indeed, ii 1 < i < m and a pair (Xi,y„_|_j), is matched then \Xi — Y„-^-j\ > d{Xi,dQ) and similarly for 
a pair {Xm+i,Yj), with 1 < j < n, \Xm+i — Yj\ > d{Yj,dQ). The term diam{Q)P\m + m! — n — n'\ 
takes care of the points oi X \jy that are not matched in the optimal matching oi X U X' and 3^ U 3^'. 
We apply the above inequality to X, y independent Poisson processes of intensity 1q ■ ly, and A"', y' , 
two independent Poisson processes of intensity Iqc • i/, independent of {X,y). Then X U X', 3^ U 3^' 
are independent Poisson processes of intensity v. Taking expectation and bounding the average of the 
difference of cardinalities in the usual way, we obtain the first inequality. 

Now, the second inequality will follow from the superadditive property of the boundary functional: 

LaQ(x,y)>Y,LopixnP,ynP). (19) 

Per 

This is proved as follows. Let {A,B,a) be an optimal triplet for LaqiX ,y): 

LaQ{X,y) - ^ - r,(,)|f + ^ qd{x,,dQr + ^ qd{Y,,dQr. 

ieA ieA" jeB" 
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li X E Q, we denote by P{x) the unique P G V that contains x. If P{Xi) = PiY^^i-^) we leave the term 
\Xi — yo.(,;)| unchanged. On the other hand if P{Xi) ^ P(Fg.(i))j from Holder's inequality, 

Eventually, we apply the inequality 

d{x,dQ) > d{x,dP{x)) 

in order to take care of the points in A''^ U B'^. Combining these inequalities and grouping the terms 
according to the cell P E V containing the points, we obtain that 

LaQ{X,y) > E f E \X,~Y,^,^\P+ qd{X,,dPr 




J2 qd{Y,,dPr+ Qd{Y„dPY 

> ^ Lap(A'nP,3^nP), 

Pev 

and we have obtained the inequality ([T^ . □ 

The next lemma will be used to reduce to uniform distributions on squares. 

Lemma 26. Assume L = Alp. Let ij., ^' be two probability measures on M.'^ with supports in Q and 
n > 0. Then 

ELi)Q{nfi) < EL9Q{n^') + 4ndiam{Q)P dTvi^J-, ^J.')■ 
Consequently, if f is a nonnegative locally integrable function on , setting a = f/vo\{Q), it holds 

ELoQinflg) < ELgginalQ) + 2ndmm{Q)P [ - a\ dx. 

Jq 

Proof. The functional Lgq satisfies a slight modification of property (TZp) : for all multisets X, y, Xi,yi,X2, 3^2 
in Q, it holds 

L9Q{XUXi,yuyi)<LaQ{XUX2,yiiy2)+ diam(Q)''(card(A'i) + card(A'2) + card(3^i) + card(J^2)). 

Indeed, we start from an optimal boundary matching of Lqq {XUX2,yuy2), we match to the boundary 
the points of {X, y) that are matched to a point in (<%2, 3^2)- There are at most card(A2) +card(3^2) such 
points. Finally we match all points of {Xi,yi) to the boundary and we obtain a suboptimal boundary 
matching of Lqq{X U Afi, 3^ U 3^i). This establishes the above inequality. The statements follow then 
from the proofs of Proposition |S] and Corollary [HI □ 

We will need an asymptotic for the boundary matching for the uniform distribution on the unit cube. 
Let Q = [0, 1]'' and denote 

LdQ{n) ELaQ(nlQ). 
Lemma 27. Assume L = Mp and < p < d/2, then 

hm — f-jr- = Pl^ 
where (3'^ > Q is a constant depending on p and d. 

Proof. Let m > 1 be an integer. We consider a dyadic partition V oi Q into m"^ cubes of size 1/m. Then, 
Lemma E5] applied to the measure nlij^ ^d(x) dx gives 

However by scale and translation invariance, for any q E V we have KLaqinl^i^^Q^^d) = m^PKLaQ{nm^'^lQ). 
It follows that 

Laqin) > m^-PLaQinm-"). 
The proof is then done as in Theorem [T5] where superadditivity here replaces subadditivity there. □ 
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6.2.2 General absolutely continuous measures 

Wc arc ready to state and prove 

Theorem 28. Assume L = Mp and < p < d/2. Let / : M'* — > R+ be an integrahle function. Then 



Proof. Assume first that the support of / is bounded. By a scahng argument, wc may assume that the 
support of / is included in Q = [0, 1]''. The proof is now similar to the one of Theorem [T71 For ^ € N, 
wc consider the partition Vt, of [0, 1]'' into 2^'^ cubes of side-length 2~^ . Let fc € N* to be chosen later. 
For P &Vk, ctp denotes the average of / over P. Applying Lemma [551 Lemma and homogeneity, we 
obtain 



2di 



t f f + EL{nf) > ELoQinf) > ^Lgpinflp) 

- (ELapinaplp) - 2ndi 2-''P f \f - ap 



= J2 (2"^^EiaQ("«p2-'=''lQ)-2nd5 2-'=f f \f - ap 



Setting as before fk = J^peVk '^p^p h{t) = LQQ{t)/t d where LgQ{t) = ELgQ{tlQ), the previous 
inequality reads as 



^ z _ 1 , p 
2n'' 



+ >^LoQ{nf) > I Mn2-'='^A.)/,!-5_2dfn5 2-^=^1 \f - f^]. 



As in the proof of Theorem 1 171 we may choose k = k{n) depending on n in such a way that lim„ k{n) = 

+0O, hm„ n^/''2~''''"' = +oo and hm„ 713 2"*^^") ( / |/— /fc(„) |) " =0. For such a choice, since hminft_i._|-oo > 
j3'^ by Lemma [27] and a.e. limj; fk ~ /, Fatou's lemma ensures that 

liminf / M-2-'=(")'^/fc(„))/,!-f > hminf / M-2-'^>)^/,(„))/,Yj > P'l I Z^"^- 

Our statement easily follows. 

Now, let us address the general case where the support is not bounded. Let t>\ and Q = [—£,£]'^. 
By Lemma [251 



2diam(Q)Mn J / + EL(n/) > ELggln/Ig). 



Also, the above argument has shown that 



n d 



We deduce that for any Q = 



Q 

Taking £ arbitrary large we obtain the claimed lower bound. □ 
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6.2.3 Dealing with the singular component 

In this section wc explain how to extend Theorem from measures with densities to general measures. 
Given a measure we consider its decomposition /i = /iac + /^s hrto an absolutely continuous part and 
a singular part. 

Our starting point is the following lemma, which can be viewed as an inverse subbadditivity property. 

Lemma 29. Let p G (0, 1] and L = Mp. Let Xi, X2,yi, 3^2 be four finite multisets included in a bounded 
set Q. Then 

LiXi,yi) < L{Xi U X2,yi U 3^2) + ^(^-2,3^2) + diam(Q)p(|A'i(Q) - 3;i(Q)| + 2\X2{Q) - 3^2(Q)|). 

Proof. Let us start with an optimal matching achieving L{Xi U A'2,3'1 U 3^2) and an optimal matching 
achieving L(A'2, 3^2)- Let us view them as bipartite graphs Gi,2 and G2 on the vertex sets (^"10^2, 3'iU3'2) 
and (^2,3^2) respectively (note that if a point appears more than once, we consider its instances as 
different graph vertices). Our goal is to build a possibly suboptimal matching of Xi and 3^i. Assume 
without loss of generality that Xi{Q) < 3^i(Q). Hence we need to build an injection from a : Xi ^ yi 
and to upper bound its cost J2x£Xi 1-*^ ~ '''(^)l^- 

To do this, let us consider the graph G obtained as the union of Gi_2 and G2 (allowing multiple edges 
when two points are neighbours in both graphs). It is clear that in G the points from Xi and 3^i have 
degree at most one, while the points from X2 and 3^2 have degree at most 2. For each x G Xi, let us 
consider its connected component C{x) in G. Because of the above degree considerations (and since no 
point is connected to itself in a bipartite graph) it is obvious that G{x) is a path. 

It could be that G{x) = {x}, in the case when a; is a leftover point in the matching corresponding to 
Gi,2. This means that x is a point in excess and there are at most \Xi{Q) + X2{Q) — (3^i(Q) + 3^2((9))| 
of them. 





Figure 1: The three possibilities for the path C{x). In blue, Gi_2, in red G2, the points in Xi U X2 are 
represented by a cross, points in 3^i U 3^2 by a circle. 

Consider now the remaining case, when C{x) is a non trivial path. Its first edge belongs to Gi,2. If 
there is a second edge, it has to be from G2 (since G1.2 as degree at most one). Repeating the argument, 
we see that the edges of the path are alternately from Gi^2 and from G2. Note also that the successive 
vertices are alternately from Xi U X2 and from 3^i U 3^2 (see Figure [2]). There are three possibilities: 

• The other end of the path is a point y £ 3^i. In this case we are done, we have associated a point 
y G yi to our point x E Xi. By the triangle inequality and since (a + b)P < + If due to the 
assumption p < \, \x — yj^ is upper bounded by the sum of the p-th powers of the length of the 
edges in C{x). 

• The other end of the path is a point y € 3^2- The last edge is from Gi_2- So necessarily, y has 
no neighbour in G2. This means that it is not matched. There are at most \X2{Q) — 3^2 (Q)| such 
points in the matching G2. 

• The other end of the path is a point x' G A'2. The last edge is from G2. So necessarily, x' has no 
neighbour in Gi^2- This means that it is not matched in G1.2. As already mentionend there are at 
most \Xi{Q) + X2{Q) - (3^1 (Q) + y2{Q))\ such points. 
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Eventually we have found a way to match the points from Xi, apart maybe \X2{Q) — 3^2(Q)| + \Xi{Q) + 
X2{Q) — {yi{Q) +3^2(Q))| of them. We match the latter points arbitrarily to (unused) points in J^i and 
upper bound the distances between matched points by diam((5). □ 

As a direct consequence, we obtain: 

Lemma 30. Let fii and fi2 be two finite measures supported in a bounded set Q. Let p G (0, 1] and 
L = Mp be the bipartite matching functional. Then 



Ei(/ii) < EL(/ii +^2) +IEL(m2) + 3diam(g)P(Vm(Q) + VJI^). 

Proof. Let A"! , A'2 , , 3^2 be four independent Poisson point processes. Assume that for i E {1,2}, Xi 
and yi have intensity measure fii. Consequently Xi U X2 and J^i U 3^2 are independent Poisson point 
processes with intensity fix + fi2- Applying the preceeding lemma [29] and taking expectations yields 

ELini) < EL{iii + PL2) + EL{n2) + 2diam(Q)P(E|A'i(g) - 3;i(g)| + E|A'2(g) - y2{Q)\)- 

As usual, we conclude using that 

E|A',(g) - y,{Q)\ < ^E{{X,{Q)-y,{Q)Y) = V2var(A',(g)) = y%dQ). 

□ 

Theorem 31. Assume that d E {1,2} and p e {0,d/2), or that d>3 and p e (0, 1]. Let L = Alp be the 
bipartite matching functional. Let ji be a finite measure on with bounded support. Let f be the density 
of the absolutely continuous part of ^. Assume that there exists a > -jz^, such that J \x\°'dfi{x) < +00. 
Then 

Moreover if f is proportional to the indicator function of a bounded set with positive Lebesgue measure 

hm =h f 



Proof. Note that in any case, p < 1 is assumed. Let us write fi — fiac + fJ-s where dfiac{x) = f{x)dx is 
the absolutely continuous part and jis is the singular part of ^. 

The argument is very simple if ji has a bounded support: apply the previous lemma with /ii = ujiac 
and /i2 — rijjLg. When n tends to infinity, observing that ^/n is negligible with respect to n^~^d , we obtain 
that 

. EL{n^iac) ^ ,. . „ EL(np) EL(ri^s) 
nm mt — 5 — < lim mt — - — p h iim sup — 5 — . 

" n^^d n n^^d n n^~d 

Observe that the latter upper limit is equal to zero thanks to Theorem [TO] applied to a purely singu- 
lar measures. Eventually liminf„ ^^("^^""i > Ld/^~5 by Theorem about absolutely continuous 
measures. 

If / is proportional to an indicator function, we simply use scale invariance and Theorem [24] in place 
of Theorem [11] 

Let us consider the general case of unbounded support. Let Q = [—£, i]'^ where ^ > is arbitrary. Let 
Xi,yi, X2,y2 be four independent Poisson point processes, such that Xi and J^i have intensity measure 
nig • pLac^ and X2 and 3^2 have intensity measure n(/is + ^Q" • Mac)- It follows that Xi U X2 and J^i U 3^2 
are independent Poisson point processes with intensity n/x. Set T := max{|z|; z £ Ai U A2 U J^i U 3^2}- 
Applying Lemma [29] gives 

^(^-1,3^1) < L{Xi U ^-2,3^1 U 3^2) + ^(^-2,3^2) + CpTP {\card{Xi) - card(3^i)| - |card(A'2) - card(3^2)|). 

Taking expectations, applying the Cauchy-Schwarz inequality twice and Lemma [?] (note that a > 2p) 
gives 

ELinflg) < EL(n/i) +Ei(n(/i, + Iqafi,,)) + Cp^J E[T^p] (^^/2nHaciQ) + i/2n(/x,(M'^) + /x,e(g^)) 
< EL{nfi) +EL{n{fi, + Igo/iac)) + c^n^+i 
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Since a > we obtain 

d—2p 

hmmf — ^^-Y^>liminf \% - limsup —f^ > 13'^ I I " ~ Pl f / " 



Q 

where we have used Theorem [^S] for the lower hmit for bounded absolutely continuous measures and 
Theorem [TH] for the upper limit. Recall that Q = [—l^Pf. It remains to let I tend to infinity. □ 

Actually, using classical duality techniques (which are specific to the bipartite matching) we can 
derive the following improvement of Lemma [30l which can be seen as an average monotonicity property: 

Lemma 32. Let p e (0, 1] and L = Mp. Let and 112 be two finite measures supported on a bounded 
subset Q C M'*. Then 



EL(/ii) < ELifii + ^2) + 3diam(g)P(VMi(Q) + VMO))- 

Proof. Since p e (0, 1], the unit cost c{x, y) := \x — y\P is a distance on M''. The Kantorovich-Rubinstein 
dual representation of the minimal matching cost (or optimal transportation cost) is particularly simple 
in this case (see e.g. [TUIIIBIIII]): for {xi, . . . ,Xn}, {yi, ■ ■ ■ ,yn} two multisets in Q, 

L({a;i,...,x„},{yi,...,y„}) = sup ^ f{xi) ~ fiy,), 

where Lipj^ q denotes the set of function / : Q — > M which are 1-Lipschitzian for the distance c{x, y) (hence 
they are p-H61derian for the Euclidean distance) and vanish at a prescribed point xq Q. Observe that 
any function in Lipj^ q is bounded by diam((5)P pointwise. 

Let X — {Xi, . . . , Xn-^} and y = {Yi, . . . , Yn^} be independent Poisson point processes with intensity 
yU of finite mass and supported on a set Q of diameter D < +00. By definition, on the event {A'^i < N2}, 

L{x,y) = inf L({x„i<j<iVi},{r„jeA}) 

Ac{l,...,N2};cai:d(A)=Ni 

inf sup V f(X,) - V f(Yi) 

> sup I ^ fix,) YI /(^^) I - ^'1^1 - 

where we have used Kantorovich-Rubinstein duality to express the optimal matching of two samples of 
the same size and used that every / g Lipi^g satisfies |/| < pointwise on Q. A similar lower bound 
is valid when A^i > 7V2. Hence, taking expectation and bounding E\Ni — A''2| in terms of the variance of 
the number of points in one process, one gets 

EL(/.) > E sup I ^ fix,) - Y /(^^) I - ^'^2R. (20) 

A similar argument also gives the following upper bound 

EL(^) < E sup I ^ fix,) - I + D'VWl (21) 

/eLipi y,<Ni j<N2 / 

Let Xi, A2, 3^1, 3^2 be four independent Poisson point processes. Assume that for i € {1, 2}, A'^ and yi 
have intensity fit. As already mentioned, Xi U X2 and U 3^2 are independent with common intensity 
A'l + /^2- Given a compact set Q containing the supports of both measures, and xq G Q we define the set 
Lipi Q. Using dini). 



Ei(/ii + fi2) = EL(A'i U ^-2,3^1 U 3^2) 

> E sup I ^ /(xi) - Y fM + E /(■^2) - Y f(y^) I - ^V2|mi + m2| 

i6Lipi,o yieyi X2GX2 y2ey2 
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Now we use the easy inequality E sup > sup E when E is the conditional expectation given Xi,yi. Since 
('^2,3^2) are independent from (A'i,3^i), we obtain 



> E sup I ^ /(xi)- J2 E E 

= E sup I ^ - ^ /(yi) 



/eLipi Q 



> EL{^_ll)-DP^/2\^I^\, 

where we have noted that the inner expectation vanishes and used (PT|). The claim easily follows. □ 
6.3 Euclidean combinatorial optimization 

Our proof for the lower bound for matchings extends to some combinatorial optimization functionals L 
defined by . In this paragraph, we explain how to adapt the above argument at the cost of ad-hoc 
assumptions on the collection of graphs (^rOneN- As motivating example, we will treat completely the 
case of the bipartite traveling salesperson tour. 

6.3.1 Boundary functional 

Let S CM."^ and e,p> 0. Set q ~ 2^^^ A 1. In what follows, p is fixed and will be omitted in most places 
where it would appear as an index. Given multisets X = {-^i, . . . , Xn} and y = {Yi, . . . , y„} included 
in R'*, we first set 



Geg, 

where 



Lis AX, y) = min \ ds,eAX^, Y,) ) , 



\x - y\P if X, y £ S, 



ds,eAx^y)-{ q(disi{x,S^)P + eP) if x&^S^yis ^^^^ 



g(dist(y, S'^Y + e^) if y £ S, x S 
Now, if X and y arc in S, we define the penalized boundary functional as 

Las AX, y) = ^min ^ L%sAX U A, U B), (23) 

where the minimum is over all multisets A and B in 5''^ such that card(A'U A) = card(3^Ui3) > kq. When 
e = we simply write Lqs- The main idea of this definition is to consider all possible configurations 
outside the set 5 but not to count the distances outside of 5* (from a metric view point, all of S'^ is 
identified to a point which is at distance e from S). 

The existence of the minimum in is due to the fact that L^Qg{X A,yvj B) can only take finitely 
many values less than any positive value (the quantities involved arc just sums of distances between points 
of X,y and of their distances to 5^). Notice that definition (P^l) is consistent with the definition of the 
boundary functional for the matching functional Mp, given by (jl8p . If X and y are independent Poisson 
point processes with intensity v supported in 5* and with finite total mass, we write LqsA^) for the 
random variable LqsA^ ,y)- Also note that ds,o,p{x,y) < \x — y\P. Consequently if card(A') = card(3^) 
then 

LlsiX,y)<LiX,y). (24) 
The next lemma will be used to reduce to uniform distributions on squares. 

Lemma 33. Assume (A1-A5). Let be two probability measures on with supports in Q and 

n > 0. Then, for some constant c depending only on k, kq, 

E,LgQ{n^) < E,LQQ{n^') + 2cndiam{QY d^yi^i, IJ.')- 
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Consequently, if f is a non-negative locally integrahle function on , setting a ^ Jq f /vo\{Q), it holds 

ELoginflQ) < ELaqinalQ) + cndiam{Q)P [ \f{x) - a\ dx. 

JQ 

Proof. The functional Lqq satisfies a slight modification of property {TZp) : for all multisets X, A'l , J^i , A2, 3^2 
in Q, it holds 

LsQiXUX^^yuyi) < i:aQ('^UA'2,3^U3^2)+Cdiam(Q)P(card(A'i)+card(A'2)+card(3;i)+card(3^2)), (25) 

with C = C'{k, kq). The above inequality is established as in the proof of Lemma [23l Indeed, by linearity 
and symmetry we should check (jisp and (jl7p for Lqq. To prove (jl5p . we consider an optimal triplet 
(G, A, B) for {X, y) and apply the merging property (A4) to G with the empty graph and m = 1 : we 
obtain a graph G" and get a triplet {G",A,B\J {b}) for {X U {a},y), where b is any point in dQ. To 
prove (fT7)) , we now consider an optimal triplet (G, A, B) for {X U {a},y) and move the point a to the a' 
in dQ in order to obtain a triplet (G, A U {a'}, S) for {X, y). 

With (|25|) at hand, the statements follow from the proofs of Proposition [8] and Corollary |9l □ 

The next lemma gives a lower bound on L in terms of its boundary functional and states an important 
superadditive property of Lgg. 

Lemma 34. Assume (A1-A5). Let v be a finite measure on and consider a partition Q = Up^-pP 
of a bounded subset ofW^. Then, if c ~ 4k(1 + kq), we have 

cyJ:y{Rd) diam(Q)P + Ei(zy) > ELogilq • z^) > ^ ELap{lp ■ v). 

Proof. We start with the first inequality. Let X ~ {Xx, . . . , Xm},y = {Yi, • • • , Yn} be multisets included 
in Q and X' = {X,„+i, . . . ,X,„+m'}, y' ~ {Yn+i, ■ ■ ■ ,^n+n'} be multisets included in Q'^. First, let us 
show that 

ci\m + m' - n ^ n'\diam{Q)P + L{X U X' ,y U y') > LaQ{X,y), (26) 

with ci = k{1 + kq). To do so, let us consider an optimal graph G for L{X U X',y U 3^'). It uses 
all the points but |m + m' — n — n'\ points in excess. We consider the subsets Xq C X and yo C y 
of points that are used in G and belong to Q. By definition there exist subsets A,BcQ'^ such that 
card(A'o U A) = card(3^o U B) and L{X U X' ,yuy') ^ L{XoU A,yoU B). By definition of the boundary 
functional and using 



LagiXo, yo) < LIq{Xo UA,yoLiB)< L{Xo U A, J^o U B) = L{X LiX',yLi y'). 
Finally, since there are at most |n + n' — m — m'| points in X uy which are not in Xq U 3^o (i-c points 



of Q not used for the optimal G), the modified {TZp\ property given by Equation yields (pS)) . We 



apply the latter inequality to X, y independent Poisson processes of intensity Ig • v, and X' , y' , two 
independent Poisson processes of intensity Iqc • i^, independent of {X,y). Then X U X', y U y' are 
independent Poisson processes of intensity i'. Taking expectation, we obtain the first inequality, with 
c = 4ci. 

We now prove the second inequality. As above, let X = {Xi, . . . , y ~ {Yi, . . . , Yn} be multisets 
included in Q. Let G € 5^ be an optimal graph for LaQ{X,y) and A = {X„i+i,--- ,Xk}, B = 
{Yn+i, ■ ■ ■ , Yfc} be optimal sets in Q^. Given this graph G and a set S, we denote by Eg the set of edges 
{i,k + j} of G such that Xi £ S and Yj G S, by Eg the set of edges {i,k + j} of G such that Xi & S and 
Yj e S"^, and by Eg the set of edges {i, fc + j} of G such that X^ g S'^ and Yj e S. Then by definition of 
the boundary functional 

LaQ{X,y) = L''aQ{XuA,yuB) 

^ \X,-Y,\P+ qd{X,,Qr+ Yl 'ld{Yj,QT- 

{i.k+j}eE°^ {i,k+j}eEl^ {i.k+j}eE^ 

Next, we bound these sums from below by considering the cells of the partition V. If x G Q, we denote 
by P{x) the imique P £V that contains x. 
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If an edge e = {i, k + j} G G is such that Xi, Yj belong to the same ceU P, we observe that e G Ep 
and we leave the quantity \Xi — Yj\P unchanged. 

If on the contrary, Xi and Yj belong to different cells, from Holder inequality, 

1^. - Y,\p > qdix„Pix,rr +qd{Y,,p{Y,rr. 

Eventually, for any boundary edge in Eq, we lower bound the contribution d{Xi, Q'^Y by d{Xi, P{XiYY 
and we do the same for Eq. Combining these inequalities and grouping the terms according to the cell 
P G V to which the points belong, 

LaQix,y) > E f E E qdix,,dPY+ E <id{yj,dPY 

Pev \{i,k+j}eE% {i,k+j}eEl, {i,k+j}eEj, 

For a given cell P, set A' ^ {X U A) n P" and B' = [y \J B) f] P" . We get 

E \X^-Y,\^+ E 'id{X,,dPY+ E <ld{Y,,dPY 

{i,k+j}i=E°p {i.k+]}eEl, {i,k+]}eEj, 

= L%c{{x n p) u A' , {y n P) u B') > Lerix n P,y n p). 

So applying these inequalities to X and y two independent Poisson point processes with intensity i/Iq 
and taking expectation, we obtain the claim. □ 

Let Q = [0, 1]'' and denote 

LdQ{n) = ELaQ(nlQ). 

Lemma 35. Assume (A1-A5). Let Q C K.'' he a cube of side-length 1. If < 2p < d, then 

hm — —jr- = I^L, 
where > Q is a constant depending on L, p and d. 

Proof. The proof is the same than the proof of Lemma [27l with Lemma l34l replacing Lemma [25l □ 

6.3.2 General absolutely continuous measures with unbounded support 

Theorem 36. Assume (A1-A5) and that < 2p < d. Let f ■.M.'^ ^ M+ be an integrable function. Then 



Proof. The proof is now formally the same than the proof of Theorem 1281 invoking Lemmas 1331 1341 and 
[35l in place of Lemmas [26l [25] and [27] respectively. □ 

Remark 6. Finding good lower bounds for a general bipartite functional L on satisfying the properties 
i'Hp), (TZp), (iSp) could be significantly more difficult. It is far from obvious to define a proper boundary 
functional Lgq at this level of generality. However if there exists a bipartite functional Lqq on R'' 
indexed on sets Q C S.'^ such that for any t > 0, ELa((Q) (nltg) = tPELQQ{nt'^lQ) and such that the 
statements of Lemmas [311 [331 [3S1 hold, then the statement of Theorem also holds for the functional 
L. Thus, the caveat of this kind of techniques lies in the good definition of a boundary functional Lqq. 

6.3.3 Dealing with the singular component. Example of the traveling salesperson problem. 

Let p G (0, 1]. We shall say that a bipartite functional L on M.'^ satisfies the inverse subadditivity property 
(Ip) if there is a constant C such that for all finite multisets A'l , J^i , ^2 , 3^2 included in a bounded set 

Q C M^, 

^(-^1,3^1) < LiXi U ^-2,3^1 U 3^2) + ^(-^2,3^2) + Cdiam(Q)f (1 + |A'i(Q) - yi(Q)| + \X2iQ) - 3^2(Q)|). 
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Although it makes sense for ah p, we have been able to check this property on examples only for p E (0,1]. 
Also we could have added a constant in front of £(^'2,3^2)- 

It is plain that the argument of Section 16.2.31 readily adapts to a functional satisfying (Ip), for 
which one already knows a general upper limit result and a limit result for absolutely continuous laws. 
It therefore provides a limit result for general laws. In the remainder of this section, we show that 
the traveling salesperson bipartite tour functional L = Tp, p € (0,1] enjoys the inverse subadditivity 
property. This allows to prove the following result: 

Theorem 37. Assume that either d E {1,2} and < 2p < d, or d> 3 and p G (0, 1] . Let L ~ Tp be the 
traveling salesperson bipartite tour functional. Let ^ be a finite measure such that for some a > j^^^ , 
J \x\°'dfj, < +00. Then, if f is a density function for the absolutely continuous part of pL, 

Moreover if f is proportional to the indicator function of a bounded set with positive Lebesgue measure 

hm ' =h f 

All we have to do is to check property (Tp). More precisely: 
Lemma 38. Assume p G (0, 1] and L = Tp. For any set , A2 , J^i , 3^2 in 0, bounded set Q 

L{Xuyi) < L{Xiux2,yiuy2) + L{X2,y2) 

+ 2diam(Q)P (1 + |card(A'i) - card(3^i)| + |card(A'2) - card(j2)|) • 

Proof. We may assume without loss of generality that card(A'i) Acard(3^i) > 2, otherwise, L{Xi,yi) ~ 
and there is nothing to prove. Consider an optimal cycle Gi,2 for L{Xi U A'2,3^1 U 3^2)- In Gi^2, rn = 
|card(A'i) + card(3^i) — card(A2) — card(3^2)| < |card(A:'i) — card(3^i)| + |card(A2) — card(3^2)| points have 
been left aside. We shall build a bipartite tour Gi on {Xl,y[), the points of {Xi,yi) that have not been 
left aside by Gi.2. 




Figure 2: In blue, the oriented cycle Gi,2, in red G2, in black G'l 2- The points in A'l U A2 are represented 
by a cross, points in U 3^2 by a circle. 

We consider an optimal cycle G2 for L{X2,y2), where (<^27 3^2) ^-^^ the points of (^2,3^2) that have 
not been left aside by Gi^2- We define ('^2';3^2 ) C {Xi^^^y'^) as the sets of points that are in G2. Since 
card(A'{) + card(A'^) = card(3^{) + card(>'^), we get card(A';) - card(3^() = -card(A'^) + card(>'^). It 
implies that the same number of points from the opposite type need to be removed in (X'^,y'^ and 
(^"2,3^2) in order to build a bipartite tour. We fix an orientation on G1.2. Assume for example that 
card(A'2) > card(3^2)i if ^ point x G Xl^Xl^ , we then remove the next point y on the oriented cycle Gi_2 
of y'^. Doing so, this defines a couple of sets {X'{ ,y'{) C {X[,y'^ of cardinality card(A'{) A card(3^{) and 

L{x[,y\)<L{x';yi). 

We define G\ 2 as the cycle on (A'" U Xl[,y'l U 3^2) obtained from G1.2 by saying that the point after 
X e {X'{ U X!l, y'l U y'i) in the oriented cycle Gi_2 is the next point y G {X'{ U X!l , y'{ U y!l) in Gi,2. By 
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construction, G'^ 2 ^ bipartite cycle. Also, since p € (0,1], we may use the triangle inequality : the 
distance between two successive points in the circuit G'l 2 is bounded by the sum of the length of the 
intermediary edges in Gi,2- We get 

L{X[' U X^', X' U 3^^') < L{Xi U ^-2, 3^1 U 3^2). 

Now consider the (multi) graph G = G[ 2 ^ G2 obtained by adding all edges of G'l 2 and G2 . This 
graph is bipartite, connected, and points in (A'", y'^) have degree 2 while those in (A'j', ) have degree 4. 
Let k be the number of edges in G, we recall that an culerian circuit in G is a sequence E = {ei, - ■ ■ , Cfc) 
of adjacent edges in G such that is also adjacent to ei and all edges of G appears exactly once in the 
sequence E. By the Euler's circuit theorem, there exists an eulerian circuit in G. Moreover, this eulcrian 
circuit can be chosen so that if = {ui^i,Ui} G G2 then e^+i = {7/^+1,?/^} G G'l 2 with the convention 
that Ck+i = ei. 

This sequence E defines an oriented circuit of points. Now we define an oriented circuit on (A'", 3^f ), 
by connecting a point x of (A"", 3^") to the next point y in (A'", 3^") visited by the oriented circuit E. 
Due to the property that G G2 implies e^+i G G, if x G X{' then y G 3^" and conversely, if a; G 3^" 
then y G X{' . Hence, this oriented circuit defines a bipartite cycle Gi in {X{',y['). 

By the triangle inequality, the distance between two successive points in the circuit Gi is bounded 
by the sum of the length of the intermediary edges in E. Since each edge of G appears exactly once in 
E, it follows that 

LiKy'i) < L{^i u X2,yi u 3^2) + ^(^-^,3^2)- 

To conclude, we merge arbitrarily to the cycle Gi the remaining points of (A'i,3^i), there are at most m 
of them (regularity {Tip) property). □ 



7 Variants and final comments 

As a conclusion, we briefly discuss variants and possible extensions of Theorem [21 For d > 2p and when 
II is the uniform distribution on the cube [0, 1]'', there exists a constant Pp{d) > such that almost surely 

lim n---Hlp{{Xi,...,X„},{Yi,...,Y„}) = Pp{d). 

n— ^00 

A natural question is to understand what happens below the critical line d = 2p, i.e. when d < 2p. For 
example for d = 2 and p = 1, a similar convergence is also expected in dimension 2 with scaling \/n\iin, 
but this is a difficult open problem. The main result in this direction goes back to Ajtai, Komlos and 
Tusnady [1]. See also the improved upper bound of Talagrand and Yukich in |17| . In dimension 1, there 
is no such stabilization to a constant. 
Recall that 



1 \ / 1 1 



71 / \ n ^ — ' n 

i=l i=l 



where Wp is the Lp-Wasserstein distance. A variant of Theorem[3]can be obtained along the same lines, 
concerning the convergence of 



n-^Wp 



where /z is the common distribution of the X^'s. Such results are of fundamental importance in statistics. 
Also note that combining the triangle inequality and Jensen inequality, it is not hard to see that 



1 (- E '5^. ' ^ lE^i (- E ' - E ^y^) ^ (- E -^x. , m) , 

i—1 i—1 i—1 



(similar inequalities hold for p > I). Hence it is clear that the behaviour of this functional is quite close 
to the one of the two-sample optimal matching. However, the extension of Theorem [2] would require 
some care in the definition of the boundary functional. 

Finally, it is worthy to note that the case of uniform distribution for L = Mp has a connection with 
stationary matchings of two independent Poisson point processes of intensity 1, see Holroyd, Pemantle, 
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Peres and Schramm [B]. Indeed, consider mutually independent random variables {Xi)i>i and iYj)j>i 
having uniform distribution on Q = [—1/2, 1/2]'*. It is well known that for any x in the interior of Q, 
the pair of point processes 



converges weakly for the topology of vague convergence to (^i, ^2), where 5i and ^2 are two independent 
Poisson point processes of intensity 1. Also, we may write 



where (t„ is an optimal matching. Now, the fact that for < p < 2d, lim„ n^^^KMp{{Xi}2^-^, {lijf^i) = 
Pp{d) implies the tightness of the sequence of matchings (t„ and it can be used to define a stationary 
matching a on (Si,S2), see the proof of Theorem 1 (iii) in [B] for the details of such an argument. In 
particular, this matching a will enjoy a local notion of minimality for the Lp-norm, as defined by Holroyd 
in [5] (for the Li-norm). See also related work of Huesmann and Sturm [7]. 
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